Creation of entangled states of distant atoms by interference 
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We propose a scheme to create distant entangled atomic states. It is based on driving two 
(or more) atoms with a weak laser pulse, so that the probability that two atoms are excited is 
negligible. If the subsequent spontaneous emission is detected, the entangled state is created. We 
have developed a model to analyze the fidelity of the resulting state as a function of the dimensions 
and location of the detector, and the motional properties of the atoms. 
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OS ' I. INTRODUCTION 



The preparation of entangled atomic states is one of the goals of Atomic Physics and Quantum Optics. These 
/■""s ' states are a key ingredient for studying some fundamental issues of Quantum Mechanics |Q, as well as for certain 
applications related to Quantum Information || . Methods proposed so far to "engineer" entanglement between atoms 
in the laboratory are based on achieving and controlling an effective interaction between the atoms that are to be 
entangled. Typically, these interactions are mediated by the electromagnetic field. For example, in cavity QED, two 
atoms can be entangled if they both interact with the same cavity mode S. This coupling of the two atoms to the 
field mode can be simultaneous or sequential (that is, one atom interacts first with the cavity mode, and then the 
other one). With trapped ions, entangled states can be produced by using the Coulomb repulsion between the ions, 
together with some laser couplings fjj. With these methods, it is always necessary that the atoms interchange some 
C particles (photons) or that they are very close to each other. 

In this paper we propose a scheme to prepare entangled atomic states using a different approach. In particular, the 
entangled state is not produced by an effective interaction between the atoms, but rather by an interference effect 
and state projection accompanying a measurement. Imagine that we have two atoms A and B, situated in distant 
' locations, both in an excited state |0). These atoms may decay to the state |1) due to spontaneous emission, producing 
one photon. A detector is placed at half the way between the atoms. After some time, if the detector clicks and we 
cannot distinguish from where the detected photon came, we will have produced an entangled state 
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where is a fixed phase. Entanglement is then achieved as a consequence of two facts: first, the impossibility to 
determine from the detection event which atom emitted the photon; second, the projection postulate in Quantum 
Mechanics, which indicates that after the detection the state of the atoms is projected onto the one which is compatible 
with the outcome of the measurement. The first effect is precisely the one that would give rise to interference fringes 
at the detector position if one would repeat several times the experiment, as it has been shown by the NIST group 
at Boulder [^[(6). The second effect has been used, for example, in the preparation of non-classical states of a cavity 
mode @. Using this method to prepare entangled states, the atoms do not need to interact, and no interchange of 
particles (photons) is required. In fact, the entanglement can be produced (in principle) in a time which is half the 
distance between the atoms divided by the speed of light. 

In practice, the method described above might not be very useful. First, it is very unlikely that the photon emitted 
by one of the atoms is detected. Second, and more important, even if one photon is detected, the second atom will 
eventually decay to the ground state thus yielding the state which is not entangled. Here we will analyze 

in some detail how an experiment can be performed in a realistic set-up. The idea is to use two three-level atoms 
with a Lambda configuration (see Fig. [I]). The states |0) and |1) are the two ground states, so that once the state 
(|l|) is prepared, it will stay. Both atoms are initially prepared in the state |0). The excitation is achieved by using 
a very short laser pulse, which (with a small probability) excites one of the two atoms to level |2). If following the 
excitation a spontaneously emitted photon is detected, an entangled state of the two atoms will be produced. The 
method presented here seems particularly timely, in view of the spectacular experimental progress reported by the 
NIST group of observation of interference fringes of the light emitted by two independent atoms [gj. In fact, the same 
experimental setup could be used to prepare atomic entangled states using our proposal. 
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In order to estimate the conditions that must be fulfilled to create an entangled state, we have developed a theoretical 
model describing the whole process of laser excitation of the two atoms, spontaneous emission of a photon, and 
detection. The idea is to represent the detector as a collection of atoms, and then to use master equation methods 
to describe the projection occurring when a detection event is recorded. In this way, the electromagnetic field does 
not appear explicitly in the formulae, making the calculations simpler. We emphasize that the model is equivalent to 
the one in which the whole state of the electromagnetic field is taken into account at all times, and the measurement 
projects its state along with the state of the atoms. This model can be easily generalized to other situations in which 
there are more atoms present, yielding entangled states of more than two atoms. 

The paper is organized as follows: in Section II we explain qualitatively the details of our proposal and discuss the 
main results, and some of the practical problems. In Section III we present the theoretical model. In Section IV we 
obtain an analytical formula for the fidelity of the final state as a function of the physical parameters involved in the 
problem. Finally, in Section IV we discuss the results and point out some possible generalizations. 

II. QUALITATIVE DESCRIPTION 

Let us consider two atoms A and B separated by a distance 2d. Each of the atoms has an internal structure which 
can be described in terms of a three-level Lambda system (see Fig. [I]). It consists of two ground levels |0) and |1), 
and an excited state |2). A photodetector is located at a distance D from the segment connecting atoms A and 
B (see Fig. fn). The detector is sensitive to photons of wavelength Ai (and/or polarization) corresponding to the 
transition |2) — * |1), which is characterized by a spontaneous emission rate Y\. It is not, however, sensitive to the 
ones corresponding to the other transition. 

Both atoms are initially prepared in the state |0). Then, they are driven by a very short laser pulse on resonance 
with the transition |2) «-> |0). As a consequence, sometimes one of the atoms (or both) will spontaneously emit a 
photon of wavelength Ai, which might be recorded at the photodetector. Most of the times, no photon will be detected 
after a waiting time t^>T\. In such a case, the atoms are pumped back to the original state |0), and the experiment 
is repeated until the detector clicks. Once this occurs, the state of both atoms will be described by a density operator 
Pa,b- The goal is to obtain a state as close as possible to the maximally entangled state (0) where cf> is a phase that 
does not change from experiment to experiment. That is, we wish to obtain a fidelity 

F= (*\pa,b\*), (2) 

close to one. 

The physical idea is that the laser pulse prepares a superposition state of the two atoms, which apart from the state 
|0)a|0)b also contains a coherent superposition of the states 1 0) ^4 1 2) ^ and 1 2)^10) b- Detection of a photon implies 
that a transition |2) — > |1) has taken place in one of the atoms, producing a photon of wavelength Ai that is detected. 
The term 10)^10)^ will thus be projected out from the atomic state, since it is incompatible with that event (the 
state |1) of one of the atoms must be present in the atomic state). Moreover, given the fact that the detector cannot 
distinguish among photons emitted by different atoms, the superposition of the states |0)a|2)_b and |2)a|0)b will be 
transformed into a superposition of the states |0)^|l)s and |1)a|0)b, i.e., it will be close to the entangled state (Q). 

In order to obtain an entangled state close to the ideal Bell state ([!]), several conditions have to be satisfied: (i) 
First, the laser pulse has to be such that the probability of exciting both atoms to the state |2)a|2)b has to be much 
smaller than the probability of exciting the relevant coherent superposition. Otherwise, it may happen that although 
we detect a photon emitted by one of the atoms, the other atom also emits a photon albeit in another direction which 
is not detected; this would spoil the fidelity F since the final state of this process would be |1)a|1)_b- In order to 
avoid this problem one must use a sufficiently weak or short laser pulse. In that case, the probability of exciting 
two atoms e 2 is of the order of the square of the probability of exciting only one atom ~ 2e. By choosing e <C 1 
one avoids the two-atom excitation. Notice, however, that the laser beam cannot be too weak since it would take 
a very long time to detect one spontaneously emitted photon, given that the detection probability is proportional 
to e. (ii) Second, the detector has to be sufficiently small. At each point of the detector the phase tfi will have a 
different value spoiling the fidelity since a detection does not specify the exact location of the event, and therefore 
the exact phase is unknown. Thus, the detector has to be such that at all points the phase is practically the same. In 
order to estimate the required size of the detector surface one can use the analogy between the situation considered 
here and the double slit experiment: the distance traveled by a photon coming from one atom or the other will be 
somewhat different at different positions, and therefore the accumulated phase depends on the position in which it 
is detected. The phase will be essentially constant over regions where the corresponding interference fringes have 
a constant visibility. Thus, the length L x of the detector along the XZ-pl&ne has to be much smaller than the 
intcrfringe distance, L x <C XiD/d. However, the detection probability is proportional to the size of the detector and 
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therefore we cannot take L x arbitrarily small, (iii) Furthermore, the dynamics of the atoms during the absorption 
emission cycle will also affect the final fidelity. In fact, every absorption or emission of photons by an atom is always 
accompanied by a recoil, which changes the atomic motional state. This leaves a trace of which atom has emitted 
the photon, thus also destroying the entanglement. In order to avoid this problem, one has to find a way "not to 
leave information about the motional states behind." This can be done, for example, by using trapped particles 
and operating in the Lamb-Dicke limit, where the recoil energy does not suffice to change the atomic motional state 
(similar to the Mosbauer effect). However, the extent to which this effect can be reduced will also depend on the 
temperature of the atoms in the trap, as well as on the propagation directions of the laser beams. 

In the following Sections we will solve in detail a theoretical model to answer all of these questions. Our result is a 
simple formula for the fidelity in which these effects are clearly separated. We consider a situation where the atoms 
are trapped in identical isotropic harmonic potentials, characterized by a frequency v and initial temperature T. We 
obtain 
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where 9\ as is the pulse area (Rabi frequency times time), and Fd yn and G geo represent a dynamical and a geometrical 
factor, respectively. More specifically, 



F goo = sine 



2XWd 2 +D 2 



where sinc(a;) = sm(x)/x. We also have 
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Here, r\ = 2-7rat p /Ai is the so-called Lamb-Dicke parameter, with at p = ^Jh/2rrw the size of the harmonic trapping 
potential ground state, T is the total spontaneous emission rate from level |2), and x 1S the angle between the 
propagation direction of the laser acting on an atom and the line that connects the atom with the center of the 
detector (we take this angle to be the same for atoms A and B). 

The first factor in Eq. (|^) accounts for the effects due to the laser excitation. That is, when 6\ &s increases, the fidelity 
decreases due to the fact that both atoms may be simultaneously excited. The geometrical factor is related to the size 
of the detector with respect to the interfringe distance. For small detectors compared with such a distance, this factor 
approaches one. Finally, the dynamical factor shows that the fidelity increases for small Lamb-Dicke parameters and 
low temperatures, and depends on the ratio v/T as well as the direction of the lasers. The highest fidelity occurs for 
cos(x) ^ 1 and rj 2 coth(7u//2fc B T) < (T/v) 2 . The first condition means that the laser direction and the direction of 
the photon emitted and recorded at the detector has to be practically the same. In that case the recoil given by the 
laser is compensated by the recoil experienced by the atom in the spontaneous emission process that is monitored at 
the photodetector, and therefore no trace of which atom has emitted is left behind. Under such circumstances avCT 
(weak confinement) is needed so that the atom does not have time to oscillate in the trap before the spontaneous 
emission takes place - this would destroy the compensation of the recoils between the absorption-emission process. 
In these limits we can approximate 
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On the other hand, under conditions of strong confinement (r <C v) although it is not possible to compensate for the 
harmful effect of the recoil by choosing the laser propagation direction, the dynamical factor can be very close to one 
in the Lamb-Dicke limit (77 -C 1). In particular, for rj 2 coth(fiv /2kBT) < 1 we have 



Fdyn — f — 2?7 coth 
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III. MODEL 



A. Master equation for the atoms and photodetector 



We consider two identical atoms A and B, centered at positions and r^, separated by a distance 2d = \tq — 
Each of the atoms has an internal structure which can be described in terms of a three-level Lambda system (see 
Fig. [|). It consists of two ground levels |0) and |1), and an excited state \2). Spontaneous emission from level |2) to 
both ground levels is possible, and is characterized by the rates To^ and wavevectors ko,i(f2), where f2 represents a 
direction and r = Tq + T\ the total decay width of the excited state. 

A detector of surface dimensions S = L X L V and efficiency tjd is situated in the Ay-plane, at a distance D from 
the segment connecting atoms A and B. The center of the detector ro and the center of atoms A and B define the 
XZ plane. Wc will describe the detector as a collection of independent point atoms located at position r, with r 
varying along the detector surface [§). These atoms have two internal discrete levels \g) and |e), which are resonant 
with the wavelength Ai = 2-k jk\. The level |e) is monitored for population at time intervals St which we will take to 
be sufficiently small so that the atomic dynamics can be neglected during that time. The level |e) has a width 7: for 
sufficiently large values of 7 our model corresponds to a broadband detector, whereas for small values it corresponds 
to a narrowband detector. The results will be independent of the specific value of 7. We will concentrate on a given 
atom C of the detector coupled to the quantized electromagnetic field, which in turn is coupled to atoms A and B. 
We will calculate the state in which those atoms are left when the atom C is found in the state |e), and we will add 
incoherently the contributions corresponding to different detection times and different positions r. In such a way we 
will be finally able to derive an expression for the density operator of atoms A and B conditioned to the observation 
of a click of the detector. 

Using standard methods of Quantum Optics, one can trace out the electromagnetic field and obtain a master 
equation for the atoms A, B and C: 



dt P 
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where C a denotes the Liouvillian superoperator describing the evolution of atom a alone, and 



i]-G (|r° - r|) {a% ® af 2 + 0% <g> of x ) p + B.C. 



J a '°P = l J ^e- ik < n >- r V2P*? ff e ik(n) - r + H.c. 



(9a) 
(9b) 



with afj — \i) a {j\ (superscripts indicate the atom, whereas subscripts indicate the states). Here and in the following 
we will use the symbol ® (tensor product) whenever we feel that it clarifies the corresponding expression. The vectors 
v A and r B are the position operators of the atoms A and B, while the vector r is treated as a c-number. The 
presence of the factor G(r) = — exp(iki |r|)/ {k\ |r|) is due to the dipole-dipole interaction (real part) and reabsorption 
(imaginary part) between atoms A, B and C, 7 giving the typical strength of this interaction. These two terms give 
rise to the excitation of atom C via a photon absorption from atoms A and/or B, which leads to a detection event. 
We have assumed k(d 2 + D 2 ) 1 / 2 3> 1, so that only the far-field part contributes to the dipole-dipole interaction. 
The Liouvillian action on atom C (detector) is given by 



C C P 



(10) 



In the absence of laser excitation, we have (a — A, B) 



C a = ±:[H? p ,p\-^(a2 2 p + pa% 2 ) 



dn 



(ii) 



+r J —N (n)e- ik ^- r \% 2 pa2 e^ n >r a 

+Tl J ^N 1 (n)e~ ik ^ r '*a? 2 pcT% 1 e lk ^ n > ra . 

Here, Ht p is the Hamiltonian describing the motion of an atom in an isotropic harmonic potential of frequency ^, and 
iVo and N\ describe the dipole emission pattern corresponding to transitions |2) — > |0) and \2) — > |1), respectively. 
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The master equation (||) can be solved formally as 

p(t) = e^ A+cB+cC ^-^p(t Q ) (12) 

J t 

This integral equation can be iterated to obtain a formal expansion in terms of S and J . Since each of these terms 
scales as l/k(d 2 + D 2 ) 1 / 2 <C 1, we can stop at the first non-vanishing order of the equation. The even terms of the 
expansion correspond to physical processes in which excitations (photons) are interchanged between atoms A and C 
(or B and C) . We have not included in Eq. (@) the (dipole-dipole) interactions between atoms A and B which would 
give rise to processes describing photon exchange, because they correspond to a very small correction of the order of 
X/kd <C 1 to the final result. Note that we should only consider the case in which atom C is detected in |e), which 
can only occur if a photon coming from A or B is absorbed; that is, the first non vanishing process in our expansion 
will correspond to the emission of a photon from atom A or B subsequently absorbed by atom C. This will give a 
contribution of the order l/k 2 (d 2 + D 2 ). Processes in which more than one photon are interchanged between atoms A 
(or B) and C, or in which (apart from the photon absorbed by C) other photons are interchanged between atoms A 
and B would give higher order contributions, at least of the order of X/k 4 (d 2 + D 2 ) 2 or l/k 4 (d 2 + D 2 )d 2 , respectively. 



B. Initial state of atoms A and B: Laser interaction 



So far, we have ignored the initial state of atoms A and B. Let us assume that they are driven by a very short laser 



pulse of duration t\ as <f 1 , v 1 . The state of atom a after the interaction is 



Pa(0) 



-ih? 



>p a (0)e 



ill? 



where p a (0) = <rg ® p£ p (Q), with 



tf p (0) (x exp(-H^/k B T) 



(13) 



(14) 



being the initial motional state corresponding to a thermal distribution at temperature T in the trapping potential, 
and e~ lhl ™ acts in the subspace span{|0) Q , |2) Q } as 



-ih? 



cos(6> las ) - isin(0ias) 



H.c. 



(15) 



Here, #i as is the rotation angle due to the laser interaction and k Q the laser wavevector acting on atom a. 

According to these equations, the effect of the laser on each of the atoms is twofold: on one hand, it excites a 
superposition of the internal states |0) and |2); on the other hand, it gives a kick to the atom. The coefficient of the 
superposition #i as can be easily varied by changing the laser intensity/duration. 



C. Detection 



We will use the following model for the detection ||. The initial state of the atom detector is \g). The evolution 
time is divided in time steps ti,t%, . . . ,t n , . . . of duration St <C T -1 , v~ l . After each time interval St, the internal state 
of atom C is measured and the state of the whole system is projected onto \g) or |e) depending on the outcome. Let 
us consider the case in which the detection at time t\,t%, ■ ■ ■ , t n has yielded the outcome \g), and the detection at 
time t n+ \ has yielded |e). To lowest order in our expansion, the unnormalized state of atoms A and B at time t — * oo 
once we have made the corresponding projections will be 

Pn = K lim e^ A+cB ^-^R(t n ), (16) 

t — >oo 

where K is a constant that only depends on 7, 7 and St, and 

R(t) = G (r A - r) af 2 p A {t)a^G (r A - r) f ® p B (t) (17) 

+G (r A r) a A p A (t) ® p B (t)a*G (v A r) f 
+ same with A <-> B, 
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with p a (t) = e cat p a {0). This expression along with other intermediate results are calculated in the Appendix. 

Since we do not know a priori at which time the detection will take place, we have to perform the sum over all the 
operators p{t n ). This sum can be transformed into an integral given the fact that 5t is smaller than any dynamical 
parameter corresponding to the evolution of atoms A and B. Moreover, we also have to integrate to all positions r 
corresponding to the detector; that is, to all positions of atom C . By doing so, we are adding incoherently all the 
contributions coming from detections at different points of the detector. Finally, we have to trace over the motional 
states of atoms A and B. The result, properly normalized, will give the averaged density operator provided the 
detector has performed a click (i.e., detected one photon). 



IV. RESULTS 



A. Density operator and fidelity 

As it is shown in the Appendix, the reduced density operator describing the internal state of atoms A and B in the 
case of detection can be written as the sum of two contributions 

p -tv(R 1 + R 2 y (18) 

Ri = cos 2 (6» las )sin 2 (6» las ) 

[M A ' A \1, 0) (1, 0| + M B > B \0, 1)(0, 1| + M A ' B \1, 0) (0, 1| + M B ' A \0, 0|] (19a) 
R 2 = sin 4 (6» las ) 
£c 
r 



where 



§M A ' A |1,0)(1,0| + ^M B ' B |0, 1)(0, 1| + ^(M A ' A + M B ' B )\1, 1| 



(19b) 



Here, we have defined 



M a,/3 = f dr ^ dtTe- rt tr tp [g (r a (t) ~ r) e lk ° p A (0) pf p (0)e^ Kr ' '<°> G (r"(t) - r) f } , (20) 

where the first integral is extended to the detector surface, the trace is taken over the motional states of both atoms, 
and ^'^(O) denote the initial motional states (|l4|). The time-dependent operators r a (t) = exp(iH" p t)r a exp(—iH" p t) 
are defined in the interaction picture with respect to the harmonic potential. 

The interpretation of Eq. ( |18| ) is very simple. The term R\ comes from processes in which only one atom is excited 
by the laser pulses and the subsequent photon emission is captured at the detector. This can be easily understood if 
one writes such a term as 

R 1 =[dr[ dtTe- Ft tr tp {m))m)\} , (21) 
Js Jo 

with 

m))A,B = G (r A (t) - r) e * kA - r »|l, 0) A)B + G (r B (t) - r) e^^^ |0, l) A . B . (22) 

The state |V , (^)) is th e superposition of two states. The first one comes from the process in which at time zero the 
laser excites atom A, including the corresponding recoil; then, at time t the atom emits a photon which is detected by 
the atomic detector at position r. The factor G (r (t) — r) includes the phase acquired during the propagation from 
the position of atom A to the detector as well as the attenuation of the probability of reaching the detector which 
is inversely proportional to the distance traveled (a solid angle factor). The second term has the same contribution 
but for the process in which atom B is excited. Since we do not take into account the exact time at which the 
photon is detected, we have to multiply \ij){t)){ijj{t)\ by the probability density that the photon is emitted at time t, 
proportional to e _rt , and integrate over time. On the other hand, since we do not know the point at the detector 
where the photon arrives, we have also to integrate the resulting expression over the detector surface, resulting in Eq. 
(pl|). Notice that retardation effects are not included in our formulation. They can be simply incorporated to this 
formula by changing t — > t — \r A ' B (i) — r|/c. Since here r A ' B and r vary over very small distances (size of the atomic 
wavepackets and detector size, respectively), the result will not be affected by retardation effects. On the other hand, 



G 



expanding the term i?2 in a similar way as Eq. ( |21| ) one can readily see that it comes from the process in which both 
atoms are excited by laser pulses, one photon emission is detected and the other not. The terms proportional to Fq 
correspond to the case in which the undetected photon is emitted in the transition |2) — > |0), whereas the ones with 
Ti correspond to the |2) — * |1) transition. 

With these expressions, we can easily calculate the fidelity (|^) as 



F= -cos 2 (0 las ) 



M B ' A t 



M A ' A + M B > B 



To . 
— su 

2r 



'las) 



(23) 



where (j> is the phase introduced in Eq. (Q). Given the fact that the size of the atom wavepackets is much smaller 
than D, we can further simplify these expressions. First, we write r a = rg + s" with |r° — r| ^> s a , the typical value 
taken by the operator s Q (of the order of the size of the atomic wavepacket). Then, we expand 



G 



(*) 



= ik Q T Q (0) 



s i(k"-ro +fei \tq —r\] 



-ikfs Q (t) ik a s a (0) 



(24) 



where k" is a vector of modulus k\ = 2tt /Ai and direction given by r — Tq . The integrals extended to the detector in 
Eq. (^p| ) can then be performed using standard methods of classical optics (substituting r by r in the denominator 
of Eq. (E4), and expanding r around r in the exponential for M A ' B and M B ' A ). Taking for simplicity \r A — r | = 



(d 2 + D 2 )V2 W e find M A,A 
M 



M B ' B = L x L y /(d 2 + D 2 ) and 
' i! - (M B < A y = M A ' A e l ^ A - r o-^ B ^) F ^ oFl 



dyn ■ 



where 



p — 



Lx/2 



d.r 



y J-L x /2 



dt< 



,-rt 



tr 



t P 



LJ2 

dye 

-Ly/2 



-ik!xd/(d 2 +D 2 ) 1/2 



(t) e -*-.s-(0) p A (0) 



tr 



(25) 

(26a) 
(26b) 



Evidently, F gco coincides with Eq. (|J). On the other hand, denoting by x the angle between k A and h A , which 
for simplicity we take to be equal to the angle between kf and k B , we obtain Eq. (||). By further choosing cf> = 
— (kyi • — ks • r'g) we obtain 



F =l cos 2 (0i as ) [1 + F gco F dyn ] + ^ sin 2 ( 



'las 



(27) 



Taking the worst case r = 0, we finally arrive at Eq. (||). 



B. Detection probability 

In order to derive an expression for the detection probability we just have to combine geometrical considerations 
with the detection efficiency t\d and the excitation probability. The probability of detection of a emitted photon is 
given by 

P ° = VD (d 2 + D 2 y/ 2 4n(d 2 + D 2 ) ' (28) 

being t/d the quantum efficiency of the photon detector. The first quotient in the expression is the cosine of the angle 
between the vector connecting the atoms and the center of the detector with a vector perpendicular to its surface. 
The second one is the solid angle extended by the detector from the atoms position. The probability that one and 
only one atom is excited and the corresponding emitted photon detected is Pq2 sin 2 ((?i as ) cos 2 (6*i as ). The probability 
that both atoms are excited and one of the emitted photons is detected is Pq2 sin 4 (6'i as ) (we neglect the process in 
which both photons go to the detector). Thus, the desired probability is 

P dct = sin 2 (g las ) ??D ^ ( ^^ 2)3/2 . (29) 
The maximum probability occurs for D = d/y/2. 
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V. DISCUSSION 



As shown in the previous sections, using our proposal, one can create states close to the maximally entangled 
state (Q). A typical test to determine whether one has succeeded or not, such as searching for violations of the 
CHSH inequalities §, would require the repetition of the experiment several times, and different measurements on 
the internal atomic states. A positive result would occur if P > 0.79, something imposing restrictive conditions on 
the parameters of the experimental setup. 

To create an entangled state of high fidelity the following conditions are required (cf. Eq. (Q) and Eq. (|^)): first, 
(i = sin 2 (6» las ) < 1; second, e 2 = dL x /[2X 1 (d 2 + D 2 ) 1 / 2 ] < 1; third, either e 3 = 2ri 2 cot\\(hv /2k B T){v /T) 2 < 1 (weak 
confinement) or £3 = 2r/ 2 coth(hv/ 2k bT) <C 1 (strong confinement). The first two conditions immediately imply a 
detection probability Pdct -C 1. In terms of these parameters we have 



2 



4 

6 



4 3 L y D ( Ai 



(30a) 

p ** = >^ 2 ZtK^) ■ (30b) 

Choosing a favorable case such as £1 = 0.1, £2 = 0.5, £3 = 0.1, still gives rise to a fidelity F > 0.8 (i.e., Bell inequalities 
are still violated). Let us analyze for this case how "distant" the atoms can be for sensible values of the parameters. 
Rewriting £2 as 



1 

62 =2 




(31) 



a value 0.5 impose D/L x ~ 50 , assuming it is not possible L x <C d ~ D that would minimize £2 while maximizing 
Pdct- Substituting in Eq. (|0tj) 

P dol = „. 8 , D i(^)\ (32 ) 

Considering an experiment is performed every 10~ 4 seconds (as it is typically the case with trapped ions) a Pdct = 10~ 4 
would correspond to a detection per second. Then, with L y /L x — 30 and a 50% efficiency, a separation of 100 
wavelengths is possible. Notice that the observation times cannot be increased arbitrarily for the deleterious effect 
caused by dark counts occurring at the detector increases consequently. 

Still, we need to asses the feasibility of £3 = 0.1 or equivalently of Pd yn = 0.9. In doing so, we will define a new 
parameter, i.e., 



^ fcl V^T' (33) 

so that rj 2 — rj 2 (T/v). The new parameter (a redefinition of the Lamb-Dickc parameter with V replacing v) allow 
us to study the behavior of Pdyn with respect to v /T . Once an atom and transition are chosen, 777 is fixed. Then, 
different values of v/T corresponds to different designs of the trap for the chosen atom and transition. In the weak 
confinement limit, for fixed 77/, £3 ~ v/T (just substitute Eq. ( |33| ) in Eq. (pb|)), whereas in the strong limit £ 3 ~ T/v. 
In both extremes, then, Fd yn approaches one. However, for the former case Eq. ( pq ) is not valid for arbitrarily low 
values of v/T unless cos(x) = 1 strictly. Any finite value of x implies Fd yn = at vjT = 0. Actually \ must be finite 
in order to avoid the laser light to impinge the detector, and therefore the best we can expect is a local maximum for 
Pdyn close to one. On the other hand, the strong confinement limit can be illusory for dipole transitions (needed to 
detection of the spontaneously emitted photon in a reasonable time). We are bounded, then, to treat Pdyn exactly. 
In Figures || the behavior of the dynamical factor with respect to v/T is displayed for two values of rji, namely, 0.05 
Figure a) and 0.3 (Figure b). The value rjj = 0.05 corresponds approximately to the case of the NIST experiment 
~§. The X angle has been set to 8°, far larger than the minimum needed to avoid the laser light to impinge on the 
detector (0.8° for D/L x ~ 50). In both figures the optimum case of sideband cooling reaching T = is compared with 
standard laser cooling at the Doppler limit and half the way to it. The maximum of v/T is set to one, corresponding 
to the trap frequency equaling a dipole transition decay rate. The value r\i = 0.3 represents in such a case a limit for 
Doppler cooling reaching £3 = 0.1. From the curves shown, Doppler cooling is far enough for guaranteeing Pdyn with 
sensible values of v/T. 
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The main problem which makes the detection probability small and prevents the creation of a macroscopic distant 
entangled state is the geometrical factor. The factor referring to the laser pulse area simply reduces by a factor of 
10 the detection probability. In order to reduce the effects of the geometrical factor, one can use lenses to collect 
photons emitted in different directions. One could also couple the atoms to optical fibers, which would allow to create 
entangled atoms over longer distances. In fact, one could embed the atoms in optical cavities, so that, with a high 
probability the emitted photons, would go to the cavity mode, and then to a fiber coupled to it. The extent to which 
this can be performed in practice depends on (near-) future developments in cavity-QED. 

One can easily generalize the scheme proposed here to the case of more atoms. For example, one can take N atoms, 
excite all of them weakly using a short laser pulse, and wait for a photodetection. In the ideal case, a state 

= -L(e* |1, 0, 0, . . . , 0) + |0, 1, 0, . . . , 0) + . . . + e^|0, 0, 0, . . . , 1)) (34) 
v iv 



would be created. By using more photodetectors and observing more detection events one could create more general 
entangled states, although with a decreasing probability of success. 
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APPENDIX: CALCULATION OF p AB 

Let us denote by S(t) the free evolution operator, i.e., 

S(t)=e^ A+cB+cC ^. (Al) 

Then, iterating twice, Eq. ( [l2| ) results in 

p(t) = S(t - t )p(t ) + [ drS(t - t) [(S a > c + J A > c )S(r)p(t ) + A «-> B] + 
Jto 

[ dr f dr'S{t - t) [{S a ^ c + J A ^ c )S{t - r') {(S A ' C + J A ' C )S(T')p(t ) + A «-> B}] + 
A^ B + 0((S a ' c ) 3 ). (A2) 

We are here interested only in its projection onto the detector atom excited state, i.e., in (e|p(t)|e). The free evolution 
of the detector atom is governed by 

e cCt a c = a c 
L "gg "gg ' 

e cCt a c = e - t7/2 cr c 

° u eg c u eg I 

e cCt a c = e- t7/2 cr c 

c u ge c u ge I 

c u ee c ee J 

and it is simply enough to be operated out of (e|p(t)|e) given the initial state p(to) = p A {to) ® P~ B (to) ® &£ g - Thus, 
(e\[S(t - t )p(t )]\e) = (e\S(t - r) [(S a > c + J a > c )S(r)p(t )} \e) = , (A3) 

so that 







(e|p(t)|e)= dr dr' e'^-^e-i^-^: 

J t J t 



(C A +C B )(t-r) 



A ^ B = 



dr / dr' e -7(*-r) e -i(r-T') x 

£q J to 



{<e| [e cBt p B (t )] ® [e^^^^^e^^'^^^^yfto)^] |e)+ 

rO^C^y^] | e) } + 



(e| [e c A (t-r) S A,c e c A r~A {to) 
A^ B. 



(A4) 



As explained in the text during the measurement process, the detector atom is projected n times onto the ground 
state before being projected onto the excited state at time Under the condition 8t <C T -1 ,^ -1 the evolution 

Liouvillians inside Eq. (A4) between t n and f n +i can be left constant so that 



dr / dr' e-^*- T) e -^- T '>i?(t„) = + e^ St - 2e^ St l 2 \ R{t n ) , 



(A5) 



where 



R(t n ) = (e\ e cBt "p B (t ) ®(e| ^^^e^ 4 "/^)^ |e) + 



(e| S^S^e^+^pftO 



\e) + A <-+ B . 



Substituting the definitions of S a C in the previous equation, changing t n by t and denoting e £ °'p Q (0) by p a (<) one 
arrives to an expression proportional to Eq. ( |l7| ) . 

To proceed further we need to integrate the free evolution of the atoms given by 



■a :fu on _a a a _a , t-i /* in — ikn(n)r a _a a _a ikn(£7)r a i n -i 

p =-i[H tp ,p \--cr 22 p - -p cr 22 +r / dile ru ' a 02 p a 20 e oy 1 + 0^1. 
In a frame rotating with the trap Liouvillian the solution results in 

p a {t) = e -i a ^ t p a (o) e -^ t + 

ct 



(A6) 



o 



dre 



-Ft 



d jj e -.k (!])r"(T) ff « pa ^ (7 « e .k (!i)r°(r) + 1 . 

Taking into the account the initial condition ( |l3| ) we have 

km e^ A+£B ^-^R(r) = [cos 2 (0 ias )a o A o ® pf p (0)+ 

9ia S K T fVe^' [ dQe-^^^^^ pg {0)e ^ B (o) e ^o(n)r- { r') + 

' d0e - ikl (fl)r B (r') e ik^r B (0) p fl (0)e -k B r B (0) ^ (Q)r fl (r') 



(A7) 



i — >-oo 

sin 2 ( 



Bin^(fl jB .)of 1 ri / dr'e 
Jo 

sm 2 (0 /as ) e - r ViG(r A (r) - r) e > kAr >y p (0)e-* kAr >)G(rV) - r) f + 
sin 2 (^.) cos 2 (tf ;as )e- r V ® ^(r^r) - r)p A (0)pf p (0)G(r B (T) - r)t + 



(A8) 



Rearranging terms, Eq. (A8) can be decomposed as Ri{t) + R2{i~) witk 

R^t) = sm 2 (6 las )cos 2 (6 las )e- rT x 

G(r A (r) -r)e lkArA(0) |l,0) + G(r B (r) -r)e- !kBrB (°)|0,l) 



PtpPtp x 
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it 



G(r A - r)e ik ' <°>|1, 0) + G(r B - r)e" ik r *(°>|0, 1) 
R 2 {t) = sin 4 (9i as )e- rr x 

{<r£ ® ^r dr' e - r ^' I ^ e -iM«)rV%* B r*(o),,z; 

G(r A (r) - r) e ikA ^(°)p t A p (0)e- ik ^^( )G(r A (r) - r)t 



ik B r B (0) »ko(n)r B (r'). 



(O)r B (r')^k B r B (0) B 



*n ® ^fi r i y d T 'e- TT ' J dfie - * 1 

G(r A (r) r)e ikAr ^°)p t A p (0) e - ik " r ^ )G(r A (r) - r)t + 
A B} . 



Pf p (0)e- 



ik B r B (0) e iki(f2)r B (r') 



(A9) 



(A10) 



Tracing over the motional states and using the cyclic property of the trace the exponential terms in i? 2 (t) cancel out 
making the integral in r' trivial. Integrating r over the detector area and r with a density re~ rr , Ri and i? 2 arc 
finally obtained. 
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FIG. 1. Sketch of the experimental setup as well as of the internal level structure of the atoms corresponding to the proposed 
experiment. 



FIG. 2. The behavior of F dyn as a function of v/F for two different values of rji and three different temperatures. To denotes 
the Doppler limit temperature, i.e., To = KF /2fcs. 
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